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ABSTRACT 


In ths present invest.igai;ion flow around a circular 
cylinder mounted on a mobile bed was studied. Two different 
tyres of sand were used to stud^?- the probable effect of grain 
size. Attempts were made to analyse the basic mechanism of 
scour, that is, horse- shoe vortex, by measuring its vorticity 
mechanically with the help of a probe. In addition to this 
data available in literature are analysed to get an insight 
in the scour mechanism* 

From the analysis of data available in literature 
it is concluded that the horse-shoe vortex gradually assumes ■ 
an elliptical shape as the scour hole develops and the 
resulting in forced 'type- of vortex. Strength and size of 
initial vortex are found to be functions of Pier Reynolds 
number. From the present vorticity measurements it is found 
that vorticity decays very fast v/ith scour development. 

Based on avai.lable data and results of the present experiments 
a vortex model is proposed for prediction of scour depth. 
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INTRODUCTION AND LITDRATURS SURVEY 

^ .1 Introduction: 

It is well known that an obstruction such as a 
bridge pier placed against the direction of flow over a 
mobile bed leads to the scouring of sediment around it. 
Knowledge of the magnitude of this scour is very important 
in safe and economic design of foundations of bridge piers. 
Till date enough data have been colleoted in laboratories 
as well as in the field v/ith a view to adequately predict the 
maximum depth of scour. However, relatively l.ittle is known 
about the mechanism of scour ; because of the processes 
involved of water and sediment movement are complicated. 

Major contributors to the understanding of the mechanics of 
scour include investigators like Tison (4), Keuthner (4), 
Posey (9), Laursen and Toch (6;, Roper at al. (11), Baker (2), 
Melville and Raudkivi (8), and Qadar (10), 

1 ,1 ,1 Description of flow field around a pier: 

The scouring process around bridge piers is due to 
horse -shoe vortex systems. The vortex filaments transverse 
to the flow in a two dimensional undisturbed velocity field, 
are concentrated by the presence of a blunt-nosed pier to 
form the horse-shoe vortex system (a blunt-nosed pier is one 
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which induces a sufficiently large pressure gradient to 
initiate the concentration of vortex filaments described 
above. All other piers are referred to as sharp-nosed). 

This concentration of vortex filaments is due to the 
pressure field induced by the pier, Figs. (1,1) and (1,2) 
have been given to aid the understanding of the horse-shoe 
vortex system. The system is not steady for all flow 
conditions studied. Schwind (1962) noted that for some 
Reynolds numbers the horse-shoe vortex is shed periodically. 
The strength of the horse-shoe vortex is determined largely 
by the geometry of the pier. To some extent the condition 
of bed also affects the strength. 

Melville (197b) found that a strong vertically 
downward flow developed ahead of the cylinder as the scour 
hole enlarged. The size and the circulation of the horse- 
shoe vortex increased rapidly and the velocity near the 
bottom of the hole decreased as the scour hole enlarged. 

1 .1 .2 Scour process; 

The dominant feature of scour process around a blunt- 
nosed pier is the horse-shoe vortex system. If the strength 
of this system is strong enough to overcome the particle^ 
resistance to motion, scour will be initiated. Sediment 
particles will be disloged free upstream of the pier and 
carried out of the scour hole either by the horse-shoe vortex 
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Fig, 1,1 * HORSE SHOE VORTEX SYSTEM 



Fig. Scour around a pier 
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system and/ or by the wake vortex system on the downstream 
side of the pier. 

As a scour hole forms beneath the vortex system, 
the vortex system sinks into the scour hole. Melville and 
Raudkivi (S) point out that the vortex strength increases 
initially as the vortex system sinl^s in the scour hole, 
because additional fluid attains downward motion. However 
Baker (2) disagrees with Melville's hypothesis of increase 
in strength of vortex. 

Both Melville and Raudkivi, and Baker observed 
the increase in size of the vortex as it sinks in the scour 
hole. As the depth of scour increases, the rotational 
velocity of the vortex system decreases, causing decrease 
in the shear stress beneath' 'the vortex system. This process 
continues till the equilibrium condition is reached. The 
equilibrium condition indicates the state when the rotational 
velocity of vortex can no longer dislodge the surface grains 
at the bed. In the case of live upstream bed, the equilibrium 
condition oscillates with the movements of upstream bed 
forms into the scour hole. 

1 .1 .3 Analyses of scouring parameters: 

Limiting our attention to the case of the isolated 
bridge pier which has been modelled as a circular cylinder, 
there are many parameters which may influence the scouring 
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phenomenon: 

Variables characterizing the fluid'- 

- g acceleration due to gravity,, 

- P density of fluid, and 

- ■y kinematic viscosity of fluid. 

Variables characterizing the bed material: 

- P density of the sediment, 

s 

- size distribution, 

“ grain form, and 

- cohesion of material. 

Variables characterizing the flow: 

“ h^ depth of approach flow, 

- U mean velocity of undisturbed flow, and 

- k the roughness of the approach flow. 

Variables characterizing the bridge pier: 
its shape , 

its dimensions, (diameter, D in case ©f circular 
pier) , 

its surface condition, and 
any protection systems, 

I Assuming that (1 ) the sediment is non-cohesive and 

\ has a uniform size, d^Q ; (2) channel is sufficiently wide 
arid with a flat bed (v;itbout dunes or ripples) ; (3) pier is 
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cylindrical j circular and perfectly smooth, the scouring 
depth h depends upon eight parameters: 

hg = f^(P, V, g, dj^^pg, h^, 0 , D) 


These parameters may he replaced by the following ones: 


hg = ^2 ^ » 13) 

with A = (p - P ) / P , the relative submerged density and 

o 

^ 1 /2 

Uc. =(gh^l)' , where I is the hydraulic gradient. 

The theorem of Vaschy-Buckinghara allows us to write: 


^s 

D V 


: {5^ 

— , 
Ag d^Q 




s 

/ * • ^ 


( 1 ) 


Sqn, (1 ) can be considerably simplified by the following 
considerations: 

- The influence of the deformation of the free 
surface on the floW' field is negligible if the Froude number 
of the flow is sufficiently low, 

- There is an emperical relation for initiation 
of motion, relating 




and 


Uc 


V As d^Q 

The term is constant by considering only 

natural sediments ( pebbles, j^ravel or sand, a '^ 1.65), 
Under these assumptions Equation (l ) may be simplified to: 
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n 



•^50 X ’^sr 
— ) or — 

D D 


h 




This means that the scour deoth h will depend 
mainly on the ratio of mean velocity to mean critical 
velocity and the relative values of j^rain size, flow depth 
and pier diameter. 

1 .2 Literature Survey: 

t.2,1 Description of model data: 

Some of the most important literature on the subject 
are summer! sed belov/: 

Chabert and Engeldinger (1956) (4) studied the 
influence of flow velocity and showed that two regimes should 
be distinguished for velocities at or below the threshold 
velocity of movement of the bed material , scour depth 
approaches a limit asymptotically (see Fig. 1.3) whereas 
for a larger velocity scour depth fluctuates due to the 
periodic dumping of material in the scour hole by moving 
dunes (Fig. 1.4). Maximum scour depth was obtained at 
vel.ocities near the threshold velocity, v/hereas scour started 
at about half the threshold velocity (Fig, 1.5). 

The influence of grain size, pier diameter and water 
depth can be seen from Fig, 1.6, which shows a small influence 
of grain size, a negligible influence of water depth for water 
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depth/pier diameter ratios larger than one, and an increase 
of scour with D in v/hich 1 and D is the diameter 

of the pier. 

Laursen and Toch (1 956,1 953)(1 6)studied the influence 
of water depth, mean flow velocity and sediment size with a 
dumb-bell pier under an angle of attack of 30^. The results 
are given in Fig, 1.7, from which it v/as concluded that 
there was no systematic influence of grain size and velocity 
in the range studied. There is an influence of v/ater depth 
as is clear from the Fig. 1,7. 

Shen, Schneider, Karaki (1966a, b, 1969) (4), H-oper, 
Schneider, Shen (1967) (4), Shen (1971) (4) • From an analysis 
of the flow field and the horse-shoe vortex system near a 
circular pier it v/as concluded that the circulation of the 
vortex is proportional to U . a ( a = b/2). The influence 
of grain size was considered to be neglegible for d ^0.5 mm. 
Roughening the upstream face of the pier to decrease the 
vertical velocities or the stren'^th of the vortex had no 
effect. 

Roper et al. (1967) (11) proposed a vortex model based 
on the assumption of potential flow approaching the cylinder 
using a control volume approach. They related the scour depth 
to Reynolds number 0 where U is the average velocity of 

the approaching flow and v is the kinematic viscosity of the 
fluid. 
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Baker, 1931 (2) proposed a vortex model based on 
the assumption that the circulation of vortex remains 
same throughout the scouring process, only the size of the 
vortex increases and correspondingly the rotational velocity 
uf the vortex decreases. Based on these considerations, 
he found that the non-dimensional scour depth varies with 
Froude's parameter U/ \l g D 

Baker's model provides a basis for semi-empirical 
correlation of the variation of the equilibrium scour depth 
upstream of a circular cylinder in a bed of sand, with the 
various flow parameters. For clear water scour 

se , \ 4 . , / h . 

= (,a^N-a 2 ) tanh ( - ) 

where , a 2 , = f n ( h/d , c } 

G ^ - 1 ) g 

r 

N =. U / [( ?-^ - 1 ) gd 

P 

Qadar (1 1 ) the rotational velocity and size of the 

vortex before scour and related these to the approaching 
flov/'s average velocity and obstruction width respectively, 

Gupta , (5) observed that the first separation point on 
the plate upstream of the cylinder occurs around r/D o; 1 . 5 
at 9=0*^, He postulated the presence of three major 
vortices on the plate at the front of cylinder. A method of 
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computation of strength of horse-shoe vortex from the 
pressure distribution on the plate was developed. 

Muzzammil (9) found that the diameter of the horse- 
shoe vortex varies from 0.2 D to 0.6 D. He related the 
maximum scour depth with the initial strength of the horse- 
shoe vorticies and the upstream vorticity parameter. The 
pressure distribution on the surface of cylinder before 
and after scour was studied, and it was found that the 
pressure remains the same above the bed level and decreases 
inside the scour hole indicating the recovery of pressure. 
He observed that the of the sediment in the scour 

hole corresponds to d^^ of the original bed sediment, 
indicating the armour effect in the scour hole, 

1,2,2 Description of field data: 

Inglis, 1949 (4) compared the total scour depth 
(sum of general scoured depth, scour due to contraction and 
local scour due to the piers) v;ith the Lacey regime depth: 

hslLaoey) = ( in M.K.S.) 

in which f = silt factor = 1.76 = median grain 

size in mm. The average value was 2,09 with a r.m.s. value 
of 12.9 percent. For design purposes generally a value of 
2.0 is used, therefore: 

Total scour depth = 0,95 (Q/f)"'^^. 



Arunachalam, 1965 (4) reanalysed the data and found 
that the correlation could be improved by leaving f out of 
correlation. The result was: 

Total scour depth = (2.09) 0.473 0.95 


Ministry of Railways , India, 1967, 1968 (4). 

1967 Report: For Railv/ay bridges in all parts of India 

for flow parallel to the pier, maximum scour depth occured 
at the nose and an average ratio of 'total scour depth' to 
Lacey depth of 1 .71 was found. 


1968 Report: (Scour around the piers of Ganga pul at Makameh), 
Average scoured depth at the nose of the piers gave a ratio 
of 1 . 75 , vsrhereas for an inclined attack scour along the 
sides averaged 2.15 ^gfLacey)* observed that during 

non-monsoon floods the value of the ratio was larger, possibly 
due to a lower silt content, according the report. The corres- 
ponding values from 13 observations were 2.94 (Lacey) 
the nose and 2.7 for scour along the sides of the 

piers, Measured scour depths during non-monsoon floods were 
slightly smaller, however, than during the largest monsoon 
flood. 


1 ,3 Present Investigation: 


From the above investigations two ma;ior approaches 


to the scour predictions are found. The first method mainly 
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aims at the prediction of maximum scour depth. In achieving 
this aim the attempt is to get a better prediction method 
by relating the various parameters to the scour mechanism, 
Contributors to this rnethodology are Tison (4), Keuthner (4), 
Posey (b), Laursen and Toch (6), Lacey (4), Inglis (4) and 
Arunachalum (4). In the second method, attempt has been 
made to understand the mechanics of the scouring process. 

Most of the studies attempted to relate the circulation of 
the horse-shoe vortex before scour with pier reynolds number, 
and other flow parameters relating to the circulation ratio. 
Strength of horse-shoe vortex is assumed cbnstant during 
the scouring process. Based on the concept of constant 
strength of vortex, a model has been developed for scour 
predictions. In the study of this process, the major 
contributors are Baker (2) and Qadar (11), and Muzzammil (l). 

In both. these method^^no attempt has been made to 

study the change in flow mechanism during the scouring 

V t c 

process, Purther|^odel developed by Baker violates the 
eraperical conviction that the strength of vortex increases 
during the scouring process as Iielvelli and Raudkivi's data 
predicts. Baker assumes that vorticity of the vortex decrease 
linearly with the increase in scour depth. There is a 
necessity to go into the details of understanding how 
vorticity decays during the scouring process, to achieve a 
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meaningful development of vortex scour model. This has 
been attempted in the present work. 



CT-IAPTER 2 


EXPERT r/ISNTAL METHODS OF DATA COLLECTION 

2.1 Details of Flume: 

Tv/o rectangular water flumes were used for the 
present investigation - one of 61 cm width and 3.0 m length 
v/hile other one was A5 cm wide and 6.0 ra long. The system 
of flow straighteners have been kept at the upstream end 
of the flume to get vortex free flow. The water has been 
supplied to the flumes from a overhead tank which get supply 
from the pumps. A calibrated rectangular weir has been 
used for measuring discharge passing through the flume . 

The test section for all experiments in both the flumes 
has been selected at 1 .0 m upstream from the outlet. A 
cylinder 5.0 cm diameter and 30 cm height has been used 
as pier. 

2.2 Details of Sediment Distribution: 

Flume A was filled v\rith sand of median size 0.50 mm 
and standard deviation 0.9Q1 mm. Flume 3 was filled with 
sand of median size 0.145 mm and standard deviation 0.0154 mm. 
The size distribution of sand is shown in Fig. 2,1, 




FIG.Z.1 SEDIMENT SIZE DISTRIBUTION 
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2,3 Details of Flow Conditions: 

Four experiments at different Froude numbers 
(one at incepient condition, two at pre-lncepient conditions 
and one at post-incepient condition) were conducted in 
'Flume A' .Every experiment v/as run for 4-5 hours till the 
scouring process attained an equilibrium condition. Table 2,1 
shows the details of flov^ conditions. 

Three experiments at different Froude numbers (one 
at incepient, one at pre-incepient and one at post-incepient 
condition) were conducted in 'Flume B' , The experiments 
were run till the scouring process attained an equilibrium 
condition. Table 2,2 shows details of flow run in 'Flume B' . 

TABLE 2.1: DETAILS OP FLOW COWDITIOMS IN FLU?® A 



48.52 


0.383 3.09 


TABLE 2.2 


DETAILS OP FLOW COIOITIONS IN FLID® B. 


Flume B 

'^50 ~ 
cTdi = 

0.145 mm 

0 . 01 5 mm 

Average Temp = 27*^8 

Kinematic viscosity, v 

=8,0x1 0“^ 
m2/sGc 

Run 

No. 

Sym- 

bol 

0. in 
m^/ sec 
x1 0"3 

cm 

6 

U 

c/s 

0 

o 

c/s 

Fr 

xIO^ 

5 

o- 

2.34 

S.3o 

6.0 

6.8 

18.53 

0.11 

1 ,36 

6 


3.28 

6.52 

5.4 

10.93 

22.15 

0.14 

1 .38 

7 


9.10 

8.47 

6.3 

23.88 

28,01 

0.24 

1 .75 

2.4 

Scour Measurements 

with 

Time: 





The scour depth has been measi red from the starting 
of scouring process to the moment v/hen it appeared to attain 
an equilibrium condition, after interval of 1 ,2 , 0,20,30, 

60,120, 180, 240, 300 minutes. Fig. 2.2A shows the plots 
between scour depth and time for the first four runs, Fig.2.2B 
shows scour depth plotted against time for the last three 
runs. The rates of scour with time (dhis/dt) at different 
time intervals , stated above, were computed from these plots 
and have been shown in Table 2.3 and Table 2,4. The values 
of dimensionless scour rate (non-dimensionalized by U^)are 
also shown in the table. 
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Time tn minutes 

riG. 2-2 SCOUR DEPTH PLOTTEO AG A INST TIME 
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2.5 Details of Velocity Measurements at the U/S: 

The velocity measurements of the flows at about 1 m 
u/ s of the cylinder alon/^ the line of symmetry were recorded 

with the help of a pitot tube. The shear velocity, Vx { = > 

where is the bed shear stress) has been computed assuming 
lof'j;arithmic velocity distribution namely, 
u ^ 

_ — 111 y constant 

v/here y is the ordinate measured from the bed, X is Karman's 
constant taken equal to 0.4. 

Figure 2.3 shows the details of shear velocity 
computations for all the runs and details of velocity 
measurements are shovm in Table 2,5. 

2.6 Velocity and Pressure Distribution in the Dcour Hole: 

Measurement of velocities within the scour 

hole at equilibrium condition, was done with the help of a 

mini-current meter along the centre line of the flume as 

well as along the depth of flow. Minicurrent meter (Type 

No. 7625), developed at Central 4ater and Power Research 

Station Poona, India, is capable of measuring the speed of 

water in the range of 5 ciVsec to 100 cn/sec. The propeller 

having 5 vanes, revolves at a speed linearly proportional to 
iir f of'-' '5 

.t-ha-t of flowing water. The revolutions are counted by/ electro 


mechanical counter. 
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Pressures within the scour hole were measured with 
the help of a probe, details of which have been shown in 
Fig. 2.4. The average value of pressures measured at a 
distance of about 1 .0 m upstream of the cylinder has been 
used as a reference pressure. 

Figs. shows the velocity contours velocities 

having been non-dirnensionalized . 

2.7 Details of Vorticity Probe, ^orticity Measurement: 

Details of the vorticitv nrobe used in the experiments 
are shovm in Fig. 2.4. The vanes are made of very thin foil 
so that their inertia is minimum possible. The friction at 
the hinges has been reduced to its minimum. The number of 
rotations per minute (r.p.m," made by vane assembly (rotor) 
have been counted by naked eyes after an interval of 1 ,5 , 

10, 20, 40, 60, 120, 180, 240 minutes. These counts have 
been directly adopted as vorticily readings. These rotations 
were recorded at the location where rotor was judged to rotate 
fastest, i.e,, centrally located in the horee-shoe vortex. 

Very small and thin pieces of flourescent tape wdre pasted on 
the same side of both the vanes so that it could be possible 
to count the number of rotations, seeing from above the water 
surface , 

The details have been given below. 
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Analysis of vorticity of vortex, as it sinks with the 
scour hole development has been analysed in the next chapter. 
Prior to this, analysis vortex strength and size before 
initiation of scour has been made from the data available in 
literature , 

2.S Sediment Size Distribution in the Deepest Portion of 
the Scour Hole*. 

The sediment grains have been collected from the deep 
portion of the scour hole by pressing the plastic clay on 
the scoured surface. Sieve analysis has been carried out 
for these sediment samples. Sieve analysis data for three 
runs on sand of median size 0,5 mm and the original bed 
sample sieve analysis data for both the sands have been 
plotted in Fig. 2.1. 



CHAPTER 3 


CHARACTERISTICS OF VORTEX 


Characteristics of vortex studied are the position 
of vortex in scour hole, vorticity decay and strength of 
vorticity before tne initiation of the scour. These 
characteristics of vortex are studied using the data available 
in literature. Data of Melville and Raudkivl, C,J. Baker, 

A.k, Gupta, Qadar were analysed. The vorticity data collected 
from tv/o median size sands beds are analysed in the light 
of the analysis of the data available in literature. Static 
pressure distributions and velocity distributions in the 
scour hole are also presented and discussed. 
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3.1 Analysis of ^ata Available in Literature: 

Characteristics of vortex, like positions of centre 
of vortex , scales of the major and minor axis of the 
sllipitieal shaped vortex, in a scour hole are computed using 
Melville and Raudkivi ' s data. 

3.1.1 Positions of centre and scales of vortex in a scour hole: 

Fig. ( 3,1 ) shows the magnitudes and directions 

of mean velocity past a 5,08 cm diameter cylinder in the 
vertical plane of symmetry ahead of the cylinder for the 
initial, intermediate and equilibrium conditions as obtained 
by Melville and Raudkivi, These data were obtained from a 45 ci 
wide glass sided flume filled with sand having a uniform 
grading curve with d^^ = 0.30 nrn, dj^Q = 0.385 mm, dg^ =0.500 mm 
and a specific gravity of 2.65, A flow of 17.12 1/ s was 
used and it yielded on the bed slope of = 0.0001 a uniform 
flow depth = 0.15 ra and a mean approach velocity 0.25 m/s. 
The shear velocity, Vh was 1.21 cm/ sec. They found that 
increases from the initial flat bed value at P = 7.91x10 
m^/s to r = 9.02x10 ^ ra^/s (at hg/D = 0.923) for the inter- 
mediate sfcour to r = 1,19x10 ^ m^/ s (at hg/0 = 1.32) for the 
equilibrium scour hole in the plane of symmetry ahead of the 
cylinder. 

From the Fig. (3.1 ), it is observed that the vortex 
plunges down into the scour hole, with its tangential velocities 
along the different elliptical paths remaining the same, as 
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scour progresses. The horse-shoe vortex assumes an elliptical 
shape in the scour hole. The major axis of this elliptical 
vortex stretches with the development of scour, whereas the 
minor axis remains fairly the same. 

The position of the centre of the vortex in the scour 
hole is measured and given in Table (3.1) . 'X’ is the 

perpendicular distance of the centre from a surface of 
cylinder,* 'Y' is its depth below the initial bed level and 
'Z' is its distance along the major axis, from the point 
where the sloping armoured bed meets the cylinder. 

TABLE 3.1 i POSITION OF CENTRE OF VORTEX 


Condition 

X 

(cm) 

X/D 

Y 

(cm) 

Y/D 

Z 

(cm) 

Z/D 

(Sm) 

hs/I 

Interme- 

diate 

3.13 

0.62 

2.34 

0.46 

3.76 

0.74 

4.68 

0.9^ 

scour 

hole 









Equilibrium 

Scour hole 

3.80 

0.74 

3.28 

0.65 

4.98 

0.98 

6.72 

1.3i 


It may be observed that the vortex centre moves away 
from the surface of the cylinder. The ratio of depth of vortex 
centre from initial bed level to the scour depth remains constai 
and is approximately 1/2. 
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By measuring the perimeter of the constant velocity contours 
of these vortices, the vorticity,(i>, has been computed as 
follows: 


For a circular vortex, the strength of vortex, 
can be written as 

P =■ wr ( 27ir) 

where r is the radius of circular path of a fluid element 
in the vortex. 


Pe 


Non-dimensionalizing and writing in terms of perimeter, 


p 2 Tf r 2 

UD 

For an elliptical vortex. 


1 "D p 2 

\ 

U ^ D 


... (3.1) 


f = (velocity of ellipical path) (perimeter of this pat 
or, n = V0 . Pg 

Substituting the value of P in the above equation (3.1) an 
expression for w can be written after simplification as 


WD 

U 


2 Tr(-^ ) 

(-£e ) 

^ D 

= 277 (!L.) 


or 


0) 



Pg/D is measured for given values of Vs/Utjjand vorticity is 
computed and tabulated as shovm in Table (3.2) below. 

TABLE 3.2: VALUES OF Oi FOR IWTER-MSDIATE AMD EQUILIBRIUM 
SCOUR HOLE 


Condition 

Velocity 

Contour 

(V/Uo) 

(1) 

Perimeter 
(in cms) 

(2) 

Velocity, U 
(in cm/ sec) 

(3) 

(3) 

w ^ — X 2 rr 

(r.p.s.' 

Intermediate 

0,2 

17,33 


1 .81 

Scour 

5.0 

Hole 

0,3 

26.67 

7.5 

1,77 


0.4 

33.33 

10.0 

1 .89 

Equilibrium 

Scour 

0*2 

34.67 

5.0 

0.906 

Hole 

0.4 

48,00 

10.0 

1 .31 


It can be concluded from the above table that an 
average the vorticity remains fairly constant along the 
different elliptical (tangential velocity) paths for a 
particular position of the vortex* This can be seen clearly 
in the case of intermediate scour hole condition. This confirms 
that the vortex is of forced type. Magnitude of vorticity 
decreases from intermediate scour hole condition to eq.uilibrium 
scour hole condition. 
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Figure (3.2) shows the tangential velocity distribution 
along major and minor axes for the intermediate scour hole 
and the equilibrium scour hole. The vsridth scale of the 
vortex is chosen when the tangential velocity starts decreasing, 
These width scales, non-dimensionalized with the pier 
diameter D, are plotted against h^/D in Fig. (3.4B). 

From the figure, it is clear that the minor axis 
scale 'b', remains fairly constant tb/D = 0.2) whereas the 
major axis scale a increases v/ith increase in scour depth. 

The increase is not linear but appears to be varying as 
square of hg/D, The following equation may be given. 


a 

D 


ao 

feP 

D 


1 

3T3 


(hg/D) 


(3.2) 


Perimeter scale of the vortex is computed and plotted 
in the Fig, (3.3), It may be observed that the growth rate 

5/2 

of vortex perimeter scale in the order of (hg/D) ' , The 

follovfing equation may be given 





hs 5/2 

D ^ 


The computed vorticity (Table 3,2) nondimensionalized 
v;ith the initial vorticity is plotted against hg/D for 
intermediate scour hole and equilibrium scour hole. The 

6 

decay of the vorticity is very fast, of the order of (hg/D) , 
as can be seen from Pig, (3.3 ). 
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Strength of vortex, non-dimensionalized with initial 
vortex strength, is plotted against hg/D in Fig, ( 3 . 4 A), It 
Increases directly with hg/D as can be seen in figure. 

Some important conclusions drawn from the ' detailed 
analysis of Melville and Raudkivi's data, have been summarised 
below. 

1 < The horse-shoe vortex sinlcs into the scour hole with 
scour progress and assumes the shape of an ellipse, 

2 , The centre of vortex is located at a half of the 
scour depth. 

3« Minor axis of the ellipse remains unchanged , while 

major axis increases very fast with the increase in 
scour depth. 

4, The vorticity along different elliptical paths of 
vortex is of forced type, 

5, Vorticity decays with increase in scour depth. 

6, Vortex strength increases with increase in scour 
depth , 

3.1 .2 Initial vortex strength computations: 

It is essential to know the initial vortex strength 
characteristics because of their importance in the scouring 
process. With this view, the data of Gupta (5), Muzzommil (9) , 
Baker (1 ) and Qadar (II) are analysed. 



F10V34B VARJAI1ON OF MAJOR MmQ*R AXE,SWH 
SCOW OE-RtN»f%- :• 
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The vortex strength is computed from the pressure 
profiles measured on the rigid bed along the line of 
symmetry upstream of the cylinder. The method followed is 
that of Gupta. It was observed that one or more depressions 
exist on the pressure profile at the location of 
vortices. This fact is utilized in the computation of 
vortex strength. As vortices are of forced types, having 
crest with edge of the vortox and through v;lth the centre 
of the vortex, strength of each vortex must bo proportional to 
the pressure difference between the crest and iliirough of a 
depression or a valley. The following formula was derived 
by Gupta (5)j 


r o 


V D U 


0 


2 Ar 
D 


1/2 

) 


( 3 . 3 ) 


where = Difference in pressure coefficients 

between edge and centre of the vortex, 

Ar = distance on the plate between edge and 
centre of the vortex, 

A c And A r have been shown in a schematic 
P 

figure below. 


t 

Cp 



V / D 
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The wind tunnel data of Gupta, for houndary layer 
depths (d/D) varyin^T from 0.5 to 2.0, are used. The data 
of Plizzamiriil ' s experiments, conducted in water flume for 
the same diameter of pier as has been used for the present 
investigation, were analysed. Baker's (l ) data, collected 
from wind tunnel experiments for laminar as well as turbulent 
flow, have also been analysed, A. Oadar ' s data collected 
from water flow in a flume were analysed. From all the data, 
fo/DU is computed and plotted in Fig, 3,5 against Pier 
'■Reynolds number where U is mean velocity of flow, D is 

pier diameter , and is kinematic viscosity. 

Here, it is necessary to check the validity of 
the computation with measured values of strength of vortex 
(iion-dimensionalized with ' TDIJ '). For this, velocity 
profiles measured by Baker and Gadar have been used. The 
computation procedure is described below, 

P'ig. (3,6) shows plots between rotational velocity 
(non-dimensionalized with maximum velocity) and non-dimensiona- 
lized length along the axis on which velocities were measured. 


Ho 


Now, — — = (Vg^) 2 n r X 


Ti DU 


7rDU 


0 D 


The values of (Vg^/U ) and (2r/D) have been shown in 
an illustrative figure drawn below. In this figure D is 
diameter of cylinder equal to 7.6 cm for Baker's data. 

D varies from 2.54 cm to 10.16 cm in case of Qadar's data. 
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FIG-3-5 relation between vortex strength 
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Fig, 3,6: Plot between rotational velocity and length of axis 
as obtained by Baker (1 ) 



50 



1 ’ 



I 

j 

1.1 



3J. 



^ Cl 



- r 

! 


— w /jay i-ia™- 

' t'-' ' 

SCI-IEMTIC DIAGRAM SHOV/ING VALUES OF (2r/D) and 
({ir/U ) 

These computed data are also plotted in the Fig. ( 3 , 5 ) 
It may he observed that the data points, computed from 
pressure profiles, concide with the data points computed 
from velocity profiles. Hence we conclude that the method 
of computation suggested by Gupta may be considered 
appropriate , 

The initial vortex strength 71 DU decreases 

with increase in Reynolds number as is cl^ar from the 
Fig. ( 3.6 ). The following empirical relations for 
DU with pier Reynolds number, may be 

suggested 

= 88.85 (U D/w) 

n DU 

This loads to an important conclusion that the 
strength of initial vortex, formed in the separating region 
before scour, is a function of pier Reynolds npraber. 

Using the velocity and pressure profiles, the size 
of the vortex is computed and plotted in (Fig.3,7) against 



S C 

aADAR's data 
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Pier Reynolds number U^D/v, It may be observed that the 
size of the vortex /evnajws co'yistavjt at 2 y/£> = 0.2 


Initial vorticity is computed from the Qadar's and 
Baker’s data as follows: 


OoD 


U 



U 


D 



where is the maximum tangential velocity of vortex and 
rQ is the radius measured from the centre. For the data of 
I^zzammil and Gupta, initial vorticity is computed from the 
strength of horse-shoe vortex before scour as 



U DU 


1 

2 


D 2 

(~ ) 

r„ 


Thus computed initial vorticity nondlmensionalized with 
D and 0 as w^D/U is plotted against Pier 

Reynolds number (VD/p) in Fig. (3,8). An empirical equation may b( 
written as 


increase 


“T™ =3.5 nxp I - J 

U 65000 


It may be observed that w^D/U decreases with 
in Pier Reynolds number of the flow. 


( 3 . 4 ) 


V/ith these studies on the vortex characteristics 
before and during the scouring process, the vorticity measuremen 
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made in the present work are analysed, 

3.2 Vorticity Decay with Scour: 

For the experiments with relatively coarse sand, 
i.e. , sand of median size 0,5 mm, vorticity decreased 
with increase in scour depth. This happened for all the 
runs of flow. But for experiments v;ith finer sand, i.e., 
sand of median size 0.1 45iMfvcrti city first increased and 
tb^n decreased in the same manner as in the case of coarse 
sand. This might have occured because of size of vorticity 
probe which was used for measuring the vorticity. In the 
initial stages, the width of the deepest portion, adjacent 
to the cylinder, for the finer sand is somewhat narrower and 
it was not possible to take the axis, of the probe well with 
in the deepest portion. In the present analysis only 
receding limb data have been considered as they correspond 
to well developed scour hole. 

Initial vorticity before the scour, could not be 
measured because once the cylinder is introduced as 
obstruction to the flow, the scouring is rapid and, hence, 
there is no time available for the measurornont of the 
vorticity, Therefore, vorticity is measured after one minute 
from the start of scour, Tfebles(Art .2 .7)summarises the measured 
vorticity data.. 

The vorticity is non-dimensionalized with pier 
diameter D and maxintim velocity as and plotted 

against hg/D as shown in Pig, ( 3,9). It can be observed 
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from the figure that the vorticity decays with increase 
in scour depth. The decay rate varies with one or more 
parameters which may he from among U^, , D, d^^, d^Q, 

^f'ter an exhaustive analysis, 
taking various possible combinations, it was found that 
the parameter , gives a better corelation. The 

D+hg 

scour rate is written as { — ) V /u^ . The scour scale 

D + hg ^ o 

fs chosen because the major axis of the vortex 

increases with scour depth from its Initial value corres- 
ponding to the size of initial vortex. The size of initial 
vortex is a function of pier diameter whereas the size of 
vortex at some intermediate scour is a function of both 

oier diameter and depth of scour. Also, to some extent 
h +D 

represents the ratio of strength of 
the vortex considered to the strength of initial vortex. 

The Threshold velocity , as per Shield's criterian, 

has been considered beoauuse it represents the initiation of 
grain motion. 

h_+D 

The vorticity wD/U^ decays with — ) Yy,c(50)'^^o 

as can be seen from Fig. ( The decay is very gradual 

h +D 

upto ( — ^ — ) \c(50)'^o ^ increases very fast with 

further increase in (“-—j — ) Y*c(50)'^^o ’ Here, it may be 

recalled that the vorticity computed from Melville and 
Raudkivi’s data also decays very fast with scour process 
as can be seen in Fig, (3.3). An empirical relationship 
of vorticity decay may bo written as follows: 




, hs . V* 

FIG.,i ^a DECAY OF VORTICITV WITH (1 + -^) - 
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WD 1 

4.32 + [ 12.88 (1 + 

Uo 

The vorticity of the initial vortex, co is found 
to he a function of \c(50)'^^o ' above equation. 

This can be computed by putting hg/D = 0 in the equation. 

The computed vorticity, non-dimensionalized with U and D 
is plotted in Fig. ( 3.8), The present study indicates 
that the vorticity values are, on an average, 0,1 times 
those obtained from Qadar, Muzzammil, Gupta and Baker's 
data for corresponding Pier Reynolds number. Initial 
vorticity data computed from Melville and Raudkivi's data is 
also plotted. This value is still lower when compared to 
the present study. This is somewhat strange and needs 
further investigation. 

5.3 Rate of Scour ? 


Scour depth against time has been plotted in 
Fig. (3.11) for all the flow runs. It may be observed that 
the scouring rate is higher in the early period; decays 
gradually with increase in time, leading to the equilibrium 
of the scour depth after about 180 minutes. 

Force equilibrium equation for the removal of sand 
grain from the scouring bed may be written as, 
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c! p r 

1 4 '■ 


(ho+O? Rd^ Pg'P^ 2 
5 ■‘^2 ^ 2 c 


nd^ P_-P dhg o 


..(3.6) 


In the above equation, the first terra represents 
the fluid force due to tangential velocity of vortex acting 




FIG.3-11 rate of scour with VORTICITY 





FIG’3-12 EQUILIBRIUM VORTICITV PLOTTED 

against maximum scour depth 



62 


In order to relate the equilibrium vorticity with 


initial vorticity u , a plot of — is made against the flow 
^^c(qc))6 ^ 

parameter in Fig, (3«13A), The' flow parameter 

''*0(0)'^ ° 

is chosen for the following reason, 

, \c(qo) chosen because of the fact that median 

size, dg^ of the sediment in the deepest portion of the scour 
hole corresponds to d^Q of the sand filled in the flume as 
observed by IVIuzzammil (q).Thc parameter \c(9)'^^^o^ 
rewritten as 




U 


Vorticity corresponding to equilibri'am condition 


Vorticity corresponding to vorticity before scour 
in the approaching flow. 


Minimum tangential velocity of the vortex corresponding 
•)<'c(90) ^ to threshold velocity of dqp, size 

.im ■ n 1 ,i,r-p 

Corresponds to major diameter of the elliptical vortex 
at equilibrium condition. 


The ratio of vorticities increases with increase 

i- \J 

in ''^^.c(qo)‘^ /^o^ approaches to unity at large values of 


V 


xc(90) 


VV- 




(Jir 


Following equation may be written as 


1 - exp [ -9,2 ] 


(3.11 ) 


3,5 Comparision of Experimental Maximum Scour Depth with 
the Computed One: 

In order to check the over all validity of equations 

derived so far in the present work, compression between 
h +D 

experimental computed val^ies of is 


being made. 
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From Eq, (3.1l), 
” = 1 - exp [ - 


9.2 !i£.iaoi'* 1 


= A , say 


or 


id D 
e 

U 


0) D 
U 


. A 


From Eq. (3.1 O), 


w D h + D p 

p /Sul t ^ 

■ I MJUa JL. *>..»UWL j ^ Q j> I Q 

U D 


(3.12) 


u D 
or e _ 

U 


0.48 

h „+D 
/ sm 


(3.13) 


equating right hand sides of Eq. 3. 12 and Eq. 5,13 


0 ) D 

A 

U 


w D 

or — = 

U 


0.48 

h +D o 
^ sm ^2 

^ D 


0.48 


A(-V") 


D 


(3.14) 


Now, from Eq, (3.4) (drived using Qadar,Mji 2 zammil, Gupta and 
Baker’s data), 

U 65000 

where R, 


1,14 

='3.5 exp [ - (•^~ ) ] 


UD 


e D 


V 


f pier Reynolds number , is vorticity 


computed from Eq. (3.4). 
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In the above equation, the value of vorticity is 10,0 


times higher (please refer Art, 3, 2) than the value of vorticity 
for the present work, hence dividing R.H.S, by 10.0, one can 
get expression for o3^D/ij, Therefore, 



0,35 Exp [ 


- ( 


R 


eP ^ \1 *1 
65000 ' 


(3.1"') 


■%jatii%g right hand side of -Iq. (3.14) and Sq. (3.15), 

0 ,48 D A A t 

' = 0.35 Exp )^*^ 


h -hD 2 1*37^ r , ^eD .1^- 

or, — ) = — Exp i. ( ■ ) 

D A 


65000 


( 5 . 16 ) 


Using above equation (h +D)/D are computed for different 

Sin 

runs and compared with experimental (h -!-D)/D in Fig. (3.1 3B). 
It may be observed that the comparision between the two is 
fiar. 

3.6 Vortex Strength Relationship 

Another model for scour prediction based on vortex 
strength at equilibrium condition was proposed by 
Gangadharaiah, l''iuzzammil and K. Subramanya (13) . In this 
model the length of major axis of the vortex at equilibrium 

hgj^-i-D 

condition is taken proportional to (“■“^- ) . Also, tangential 

velocity of the vortex should be equal to threshold velocity 
[V^q^qq^j of the sediment at the deepest portion of scour. 
Therefore, the strength of the vortex can be written as 

a ■ 5 'Vo(90) 


(3.17) 
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FIG.343CA) DECAY OF VORTICtTY 



13(B) COMPARISON OF EXPERIMENTAL AND COMPUTED 
RESULTS 
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Non-dimenslonalizing with 

Vy.- 6 6 

The vortex strength is plotted against the flow parameter 

UqD 

"" T Fig. (3.14) and the variation expressed as 
V-x o 

^ 0.84 

6 V*o ^ 


(3.18) 


Experimental data of the present work are incorporated 
in the above model as shown in Fig. 3.14. It may be observed 
that the present data deviates from the curve proposed. 


3.7 Velocity Distribution in Scour Hole Region: 

Velocity contours, in cm/ sec, are drawn in Fig. (3. "I 5) 
and Fig. (3.16) for Fr = 0.11 and 0,17 for the case of sand 
with median size 0.5 mm. The contours gradually rise to the 
top surface indicating the conversion of velocity head into 
pressure head. Velocity contours near the scour hole dip 
down into the scour hole. Details of velocity distribution 
in the deepest portion and very near to the armoured bed could 
not be recorded because of the large size of the current meter. 


3,8 Static Pressure Distribution in the Scour Hole: 

Pressure coefficient contours have been plotted in 
Figs. (3.15, 3.16 and 3.17). Static pressure on the line of 
symmetry gradually recovers as the flow approaches the cylinder. 
The contours dive down into scour hole indicating the sinking 
of vortex in the scour hole. 



Notation 














0-0 
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CHAPTER 4 

CONCLUSIONS AND RECOMMENDATIONS 
4.1 Conclusions 

The following? conclusions may be drawn from the 
present work. 

(a) Combining various inferences from the analysis of 
Melville and Raudkivi's (0) data one may conclude 
that vortex in the scour hole is of elliptical shape 
and forced vortex type, having its minor axis scale 
unaltered and major axis scale increasing with 
increase in scour depth, Vorticity along the different 
elliptical paths for a particular position of the vortex 
remains fairly constant, but it decays very fast 
with increase in scour depth. 

(b^ Analysis of the data of Baker (2), Gupta (5), Muzzammil 

(9) and Qadar (11) showj/ that the initial vortex 

strength decreases with pier Reynolds number. The 

size of the vortex decreases gradually with increase 

in pier Reynolds number and remains fairly constant 
/ 2r 

— - 0.2) at higher pier Reynolds numbers. 

(c) From the present experimental study, vorticity is found 

to be directly related to rate of scour hole development, 
which decreases with Increase in scour depth. Vorticity 
at equilibrium scour condition is a function of maximum 
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scour depth and threshold velocity of the armour coat. 

A model for scour depth prediction is developed based 
on the relationship of initial vorticity to equilibrium 
vorticity, in terms of Pier Reynolds number. 

(d) Contours of Velocity and pressure distribution in 
a region near the scour hole have been 

studied. Static pressure is found to recover as 
the flow approaches the cylinder in the scour hole 
region and also along the depth of the cylinder. 

Velocity contours are inclined into the scour hole 
indicating the sinking of the vortex in the scour 
hole , 

4.2 Recommendations' 

In the present work vorticity was measured with the 
help of a probe which gave vorticity values 0.1 times those 
obtained from the data of Qadar, Kuzzammil, Gupta and Baker 
for corresponding pier Reynolds number. This has to be 
investigated. One probable cause may be friction at the 
bearings. Use of better bearings or some still better device 
may eliminate this difference. 



73 


REFERENCES 

1, Baker, C.J, (1979), Laminar Horse-shoe Vortex, Journal of 

Fluid Mechanics, Vol,95, Part 2, pp. 347“367. 

2, Baker, C.J, (1980), Theoretical Approach to Prediction 

of Local Scour Around Bridge Piers, Journal of 
Hydraulic Research. Vol.18, No,1, pp. 1-12. 

3» Baker, C.J, (1984), The Position of Points of Maximum and 
Minimum Shear Stress Upstream of Cylinders Mounted 
Normal to Flat Plates. Journal of Wind Engineering 
and Industrial Aerodynamics, 18(1985) 263-"274, 

4. Breusers, H.N.C,, G, Nicollet, and H.W. Shen, (1977), 

Local Scour Around Cylindrical Piers, Journal of 
Hydpaulic Research, Vol.l5, No. 3, pp. 211-252. 

5. Gupta, A.K. (1984): 'Experimental Investigation: Boundary 

Layer Flow Past a Circular Cylinder Mounted on a 
Flat Plate’, M.Tech, Thesis, Department of Civil 
Engineering, Indian Institute of T echnology, Kanpur . 

6. Laursen, E.M. and Tocb, A. (1956), Scour Around Bridge 

Piers and Abutments, Proc. lAHR, Int, Hydr. Conf. 
Minnesota, pp. 123-131. 

7. Melville, B.W, (1975): 'Local Scour at Bridge Sites', 

Univ. of Auckland, School of Sngg. Auckland, Newzeeland, 
Rep, No. 117. 



74 


8. Melville, and Flaudkivi,A, J, (1977), Flow Characteristics 

in Local Scour at Bridge Piers, Journal of Hydraulic. 

9. f'^uzzammil , M, (1985): 'Experimental Investigation: 

Open Channel Flov; Past a Circular Cylinder ^founted 
on a Rigid Bed and on a Mobile Bed', M.Tech. Thesis, 
Department of Civil Engineering Indian Institute of 
of Technology, Kanpur. 

10, Posey, C.J, (1974), Tests of Scour Mechanism at Bridge 
Piers, Proc. ASCE, Vol. 10, HY 12, pp. 1773-1783. 

Qadar, A, (1981), The Vortex Scour Pfechanism at Bridge 
Piers, Proc, Inst. Civil Engrs, Part 2 , Vol, 71, 
pp. 739-757. 

Roper, A.T., Schneider, V.R, and Shen, H.W, (1967), 

Analytical Approach to Local Scour, Proc, 12th lAHR, 

Congress, Fort Collins, 3, pp. 151-161, 

Gangadharaiah, T., Muzzammil, M, and >Jubramanya, K. , 

Vortex Strength Approach for Bridge Pier Scour 
Predictions, Second International Workshop on, 

Alluvial River Problems, University of Roorkee, 

October 24-26, 1985. 


11 . 


12 . 


13. 



0>^^' CJiA 


